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INTRODUCTION 


A  formulation  of  the  eigenvalue  problem,  here  referred  to  as  the  canonical  eigenvalue  method  or 
canonical  eigenvalue  theory,  was  developed  to  facilitate  analysis  of  mode  coupling  effects  in  double 
ducts  in  underwater  acoustics.  This  canonical  eigenvalue  method  can  also  be  used  if  eigenvalue  curves 
as  a  function  of  some  variable  parameter,  such  as  frequency,  are  required  to  analyze  some  acoustic 
phenomenon.  Also,  if  one  or  more  eigenvalues  are  required  at  many  values  of  such  a  parameter,  as 
for  instance  in  broadband  simulations,  the  canonical  eigenvalue  method  is  a  useful  eigenvalue  compu¬ 
tation  technique. 

This  report  is  one  of  a  series  of  reports  that  explain  the  canonical  eigenvalue  method  and  demon¬ 
strate  its  application  to  various  acoustic  ducts,  including  some  double  ducts. 


The  chief  objective  of  this  report  is  to  present  the  general  theory  of  canonical  eigenvalues.  In  the 
approach  taken  here,  the  eigenvalue  equation  is  couched  in  purely  mathematical  terms,  involving  two 
dimensionless  variables  and  a  set  of  dimensionless  parameters  related  to  the  sound-speed  profile.  The 
physical  quantities  of  interest,  such  as  frequency  and  phase  and  group  velocity,  are  generated  after  the 
eigenvalues  have  been  determined.  This  approach  has  a  nurpber  of  advantages  over  customary  ap- 
proaches^as  will  be  discussed  in  detail. 


(<  pJ  c- 


The  basic  concept  of  canonical  eigenvalues  was  introduced  in  Ref.  1,  where  the  eigenvalue  equa¬ 
tion  for  three  simple  ducts  was  expressed  in  terms  of  one  mathematical  variable.  For  an  unbounded 
refractive  duct,  the  equation  also  involves  one  dimensionless  profile  parameter.  The  extension  of  the 
canonical  eigenvalue  approach  to  more  complicated  profiles  was  suggested  and  the  first  steps  of  the 
process  were  outlined  in  Ref.  1,  which  also  contains  a  brief  history  of  how  the  canonical  approach 
came  about. 


Section  1  extends  the  basic  concept  of  Ref.  1  into  a  definitive  treatment  of  the  general  theory  of 
canonical  eigenvalues.  The  theory  is  developed  so  that  it  may  be  applied  to  propagation  in  any  multi¬ 
layered  profile  in  which  the  squared  index  of  refraction  is  linear  in  each  layer  and  is  piecewise  con¬ 
tinuous.  The  method  is  extended  to  include  improvements  in  the  treatment  of  group  velocity  and  the 
determination  of  normalization  coefficients  for  the  eigenvalues.  Further  generalizations  of  the  method 
are  outlined,  as  well  as  its  application  to  complex  eigenvalues  or  complex  parameters  in  the  sound- 
speed  profile. 

Section  2  presents  an  outline  of  the  steps  necessary  to  evaluate  acoustic  fields  by  the  canonical 
approach.  Ref.  2  is  next  in  this  series  of  reports  on  double-duct  propagation  and  associated  mode  the¬ 
ory  and  applies  the  theory  of  canonical  eigenvalues  to  the  class  of  one-layer  bounded  profiles,  as  dis¬ 
cussed  in  detail  at  the  end  of  SerJo.  2. 
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SECTION  1.  CANONICAL  EIGENVALUE  APPROACH 

This  section  presents  the  general  theory  of  the  canonical  eigenvalue  approach. 

BASIC  EQUATIONS 

The  basic  canonical  eigenvalue  equations  are  reviewed  in  this  section.  The  derivation  of  these 
equations,  together  with  additional  explanatory  material,  may  be  found  in  Ref.  1. 

The  sound-speed  profile  is  layered  in  depth  and  there  is  no  dependence  on  range.  The  sound 
speed  in  each  layer  of  the  profile  is  expressed  as 

[C,/C(Z)  ]2=  l-2y,(Z-  Zi)/Ci  (1) 

where  Ct,  Z;_  and  y/  are  the  sound  speed,  depth,  and  sound-speed  gradient,  respectively,  at  the  top 
of  layer  /.  The  canonical  eigenvalue  method  is  limited  here  to  the  case  of  continuous  sound  speed  at 
layer  interfaces;  i.  e.,  the  sound  speed  at  the  top  of  layer  i  and  at  the  bottom  of  layer  i-1  is  the  same, 
with  a  value  of  C, . 

The  unnormalized  mode  functions  are  given  by 

F,  (Z)  =  Di  Ai(-  £;)  +  Et  Bi(-  &)  .  (2) 

For  the  profile  of  Eq.  1,  F,-(Z)  satisfies  Airy's  differential  equation.  Thus,  the  solution  is  a  linear  com¬ 
bination  of  the  Airy  functions  Ai  and  Bi,  with  Di  and  Ei  being  coefficients  independent  of  depth. 

The  argument  of  the  Airy  functions  is  given  by 

ti(Z)  =  (of(Z-Z,)  +  wl(C?-C?]/a]  (3) 

where  u  =  2ir/,  Cp  is  the  mode  phase  velocity,  and  aj  is  defined  by 

of  =  -  2 Yi  <u2/Cf  .  (4) 

The  sound-speed  gradient  at  the  bottom  of  layer  i  is  given  by 

Yio  •  Cl ,  y,/Cf  .  (5) 

Equation  3  may  be  expressed  in  terms  of  C  rather  than  Z  by 

&( Q  =  (1  -  (C/Cp)2]  (C,/C)2  |y<|*2/J  n2'1  f'2  .  (6) 

In  order  to  satisfy  the  boundary  and  interface  conditions,  the  eigenvalue  equation  will  involve  Airy 
functions  with  arguments  which  are  the  negative  of  Eq.  6  as  evaluated  at  boundaries  and  interfaces. 
Equation  6  leads  to 


til  ■  &(Z<)  =  ii(Ci)  =  It-iCi/Cp)2]*2'2/1'3  \Yi\-2'3  (7) 

for  the  upper  interface  and  to 

bo  -  b(Zu i)  =  6(C4>1)  =  [l-fCi./g1]  n2'3^3  |yi0|-2/3  (8) 

for  the  lower  interface. 
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We  next  define  a  dimensionless  profile  parameter  p,  as  given  by 


Pi  =  (-Wy/+i)1/3  •  (5) 

It  follows  then  that 

ii*  1.1  =  p?  tio  (10) 

and 

aui  =  -P ]l  at  .  (11) 

The  introduction  of  p,  is  a  customary  procedure  in  most  normal-mode  approaches  to  the  profile  of 
Eq.  1  because  Eq.  10  and  11  simplify  the  eigenvalue  equation. 

All  the  equations  thus  far  can  be  applied  to  the  customary  approach  to  generating  eigenvalues. 
This  customary  approach  consists  of  solving  the  eigenvalue  equation  to  determine  phase  velocity  as  a 
function  of  frequency  for  some  given  set  of  profile  parameters.  For  some  fixed  desired  frequency,  one 
evaluates  the  and  £,o  of  Eq.  7  and  8,  using  a  first  estimate  of  phase  velocity  Cp-  One  then  evalu¬ 
ates  the  eigenvalue  equation  for  these  values  of  £a  and  5,0  and  iterates  on  Cp  until  the  eigenvalue 
equation  is  solved.  This  iteration  is  repeated  for  each  mode.  Then  the  entire  process  is  repeated  for 
each  desired  frequency  to  eventually  yield  the  phase  velocity  vs.  frequency  curves  for  each  mode. 

At  this  point,  we  take  a  different  tack.  We  introduce  the  dimensionless  mathematical  variables  x 
and  y  as  defined  by 


*  *  (12) 

and 

y  =  tto  •  (13) 

Here  x  is  the  negative  of  the  Airy  function  argument  as  evaluated  at  the  top  of  layer  1,  i.e.,  at 
the  upper  boundary  of  the  duct;  whereas  y  is  the  negative  of  the  Airy  function  argument  as  evaluated 
at  the  bottom  of  layer  1.  As  we  shall  demonstrate  in  the  next  subsection  of  Section  1,  the  arguments 
of  all  the  Airy  functions  in  the  eigenvalue  equation,  i.e.,  -£/i  and  -  £/o ,  can  be  written  as  linear  func¬ 
tions  of  x  and  y  with  dimensionless  coefficients  that  are  independent  of  /  and  Cp.  The  canonical 
eigenvalues  are  the  curves  of  y  vs  x  which  satisfy  the  eigenvalue  equation  for  a  fixed  set  of  these 
dimensionless  coefficients. 

Assume  now  that  we  have  the  x,  y  solutions  to  the  eigenvalue  equation.  The  next  task  is  to  associ¬ 
ate  these  mathematical  variables  with  the  physical  variables/  and  Cp •  We  first  derive  an  expression 
for  the  frequency.  Consider  Eq.  7  and  8  as  evaluated  for  /'  *  1;  i.e.,  the  left  sides  are  x  and  y, 
respectively.  We  eliminate  Cp  from  these  two  expressions  and  solve  for  /  to  obtain 

/  =  (y-xmcla1- ty3'1 1 n\*-'  •  (14) 

The  requirement  that  the  frequency  be  real  is  an  aid  in  the  interpretation  of  the  canonical  eigenval¬ 
ues.  This  requirement  is 


y  >  x  -*  C\>  Ci 


(15) 
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and 


y  <x  -*•  Ci  <C2  .  (16) 

We  shall  refer  to  a  plot  of  \y  -x\3/2  vs.  r  as  a  canonical  frequency  plot,  which  is  independent  of 
profile  parameters.  The  frequency  may  be  obtained  from  such  a  plot  by  the  scale  factor 
| CjCi2  -  ir3/2|yi|jr‘\  as  evaluated  for  the  desired  parameters  of  the  first  layer  of  the  profile.  Once 
the  frequency  of  Eq.  14  has  been  evaluated,  the  phase  velocity  may  be  evaluated  as  a  function  of  fre 
quency  from  the  expression 


Cp  =  c,[i-.rJ/Jff-J/,|yilI/Jx]-|/2 .  (17) 

Equation  17  is  obtained  by  solving  Eq.  7  for  Cp  for  /  =  1.  At  first  sight,  it  would  appear  that  /  and 
Cp  only  depend  on  the  profile  parameters  of  the  first  layer.  We  note,  however,  that  /  or  Cp  depend 
on  x  and  y  or  x,  respectively.  These  variables  in  turn  depend  on  all  of  the  profile  parameters  as  con¬ 
tained  in  tlw  eigenvalue  equation.  Thus,  /  and  Cp  are  explicit  functions  of  the  first  layer  parameters 
and  are  implicit  functions  of  the  other  profile  parameters. 

FORMULATION  FOR  ARBITRARY  MULTILAYER  PROFILES 

The  major  purpose  of  this  subsection  is  to  demonstrate  how  the  canonical  eigenvalue  equation  can 
be  formulated  for  an  arbitrary  piecewise  continuous  profile  of  the  form  of  Eq.  1.  The  only  problem  to 
be  solved  is  to  obtain  a  suitable  form  for  the  and  £io  of  Eq.  7  and  8.  These  are  the  negatives  of 
the  various  Airy  function  arguments  that  appear  in  the  eigenvalue  equation.  Other  features  of  the 
eigenvalue  equation,  such  as  the  particular  boundary  conditions  and  the  values  of  p,  of  Eq.  9  that 
appear  as  coefficients  of  the  Airy  functions,  pose  no  problem  in  the  formulation. 

Our  goal  is  to  express  and  £io  as  linear  functions  of  x  and  y  with  coefficients  that  are  inde¬ 
pendent  of  /  and  Cp.  Let  us  assume  that  we  have  expressions  for  t/i  and  £/o  and  wish  to  derive  the 
expression  for  f,+I>I  and  0  in  terms  of  and  f,0.  The  solution  for  f/+1  ,  is  given  by  Eq.  10, 
which  is  used  in  all  eigenvalue  approaches  that  introduce  the  dimensionless  parameter  p,- .  Our  task  is 
to  determine  the  coefficients  Af/+1  and  Nui  such  that 

C/»i,o  =  Cii  +  Nui  &o  •  (18) 

Substitutions  from  Eq.  7  and  8  lead  to 

[1  -  (Cuz/Cp)2]  x2'3  f  2'3\yux.o\'2/)  =  (19) 

Muz  [  1  -  «VC,)  V/3  / 2/3M'2/3  +  Nux  1 1  -  (Cui/Cp)2]*2'3  f  3/3|y,0|-2/3  . 

We  see  immediately  that  the  form  of  Eq.  18  must  be  homogeneous  in  Cn  and  £,•<>.  for  the  addition  of 
a  constant  term  to  Eq.  18  would  violate  the  requirement  of  independence  of  frequency.  Any  constant 
term  would  need  to  contain  f 2/J  as  a  factor. 

For  the  homogeneous  equation,  jt2/3  f  2/3  cancels  out  of  Eq.  19  as  a  common  factor.  Since  A/,,] 
and  Nux  are  to  be  independent  of  Cp,  Eq.  19  leads  to 

|y/*i.ol'2/3  =  MiM'm  +  tf<*.|y/o|-2/3  (20) 

and 

CUyux. o|-2/3  =  MUlCj\Yl\'2/3  +  Nhl  C?Jy,o|-J'3  .  (21) 
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Equation  20  is  obtained  by  equating  terms  in  Eq.  19  that  do  not  contain  Cp ■  Equation  21  equates 
these  terms  of  Eq.  19  that  contain  (1  /Cp)1  as  a  common  factor  that  has  now  been  canceled. 

From  Eq.  S,  we  may  obtain 


|yfol-2/3  =  |y/|-2/3  cj/cu 

(22) 

From  Eq.  9,  we  obtain 

|yuil'2/J  =  p/|y»ol_2/J  • 

(23) 

From  Eq.  22  and  23,  we  obtain 

|y«M.ol’2/3  =  Pf|y/|-2/3cf/cf+2  • 

(24) 

We  next  substitute  Eq.  24  and  Eq.  22  into  Eq.  20  and  cancel  common  factor  |yj|'2/3  to  obtain 


P?  C*/Cl2  =  a/m  +  Mm  Cf/cl,  .  (25) 

We  substitute  Eq.  24  and  Eq.  22  into  Eq.  21  and  cancel  common  factors  C2|yi|-2/3  to  obtain 


Pi  =  A/j+i 

+  M*i  • 

(26) 

Thus 

A/m  =  pf 

-  Mm  • 

(27) 

We  next  substitute  Eq.  27  into  Eq. 

25  to  obtain 

Mm 

*  P?[(Q/Cm)2- 

1]/[(C,/Cm)2-  1]  • 

(28) 

We  note  that  A/m  and  Mm  do  not  depend  upon  /  or  Cp  and  depend  only  on  functions  of  dimen 
sionless  ratios  of  sound  speed  and  gradients. 

The  analysis  is  simplified  if  we  introduce  the  parameter  5m  as  given  by 


Sm  =  l( Ci/Cui)1  ~  1]/[(C//Cm)2  -  1]  •  (29) 

It  follows  from  Eq.  28  and  27  that 

Mm  =  P?  5m  (30) 

and 

A/m  =  p?(1-Sm)  •  (31) 

Equations  10  and  18  allow  us  to  determine  all  values  of  Jn  and  f/0  in  terms  of  x  and  y  by  a 
recursion  process.  We  introduce  the  coefficients  P(i,  P/0,  Qi\,  and  Qj0  as  defined  by 


and 


The  coefficients  are  given  by 


£u  =  P ii  x  +  Qi i  y 

(32) 

Cio  =  Piox  +  G<0y6. 

(33) 
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P\\  *  i  ;  Qu  =  0 

(34) 

Pio  -  0  ;  Q\o  =  1 

(35) 

^21  =  0  ;  Qi\  =  Pi 

(36) 

P 20  =  p?(l  —  ^2)  ;  Q10  -  pf  St 

(37) 

P 31  =  pipiO  -S2)  ;  Qj\  =  P2P1  Si 

(38) 

and 

P 30  =  plp?S3(l-Si)  ;  Qi 0  =  pip?  [1-Sj(1-52)]  . 

(39) 

Equations  34  and  35  follow  from  Eq.  12  and  13,  respectively.  Equation  36  follows  from  Eq.  10. 
Equation  37  follows  from  Eq.  18,  30,  and  31.  Equation  38  follows  from  Eq.  10  and  37.  Equation  39 
follows  from  Eq.  18,  36,  and  37.  This  process  can  be  continued  for  as  many  layers  as  desired. 

This  procedure  assumes  that  the  profile  has  both  upper  and  lower  boundaries.  The  same  scheme 
works  when  the  deepest  layer  is  an  infinite  half  space.  When  the  deepest  layer  is  bounded,  one 
requires  both  £(1  and  f,o-  However,  when  the  deepest  layer  is  unbounded,  only  £;i  is  required. 

The  procedure  must  be  modified  slightly  when  the  upper  layer  is  an  infinite  half  space.  Here  we 
associate  x  and  y  with  upper  and  lower  arguments  of  the  second  layer.  Thus 

^21  =  l  ;  Qn  =  0 

(40) 

and 

P 20  =  0  ;  Qio  =  1  ■ 

(41) 

Higher  order  terms  may  be  obtained  from  Eq.  36  to  39  by  increasing  the  subscripts  by  1.  It  follows 
also  that 

P\o  =  Pi*  and  Qio  =  0  . 

(42) 

This  example  points  out  that  x  and  y  do  not  have  to  be  associated  with  layer  1 
associate  x  and  y  with  any  layer,  e.g. 

.  Indeed,  we  may 

x  =  5ii 

(43) 

and 

y  =  Sio  • 

(44) 

We  may  obtain  higher  order  terms  by  Eq.  34  to  39  by  increasing  the  subscripts  by  /- 1.  We  may  obtain 
lower  order  terms  by  running  the  recursion  process  down  rather  than  up.  Equation  10  and  18,  respec¬ 
tively,  may  be  manipulated  to  yield 

S/-1.0  =  pm 

(45) 

and 
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(46) 


&-i.i  =  Mj'iZio  ~  Ni  tt-u o)  • 

We  illustrate  the  use  of  Eq.  43  and  44  with  layer  1  having  an  upper  bound  rather  than  the  infi¬ 
nite  half  space  of  our  previous  example.  Thus,  Eq.  40  and  41  hold.  Equation  45  leads  to  Eq.  42, 
while  Eq.  46  leads  to 

Pn  =  -pI2(l-S2)-‘  Si-.  Qu  =  p,'2(l-52)-1  .  (47) 

If  one  defines  x  and  y  by  Eq.  43  and  44,  C t.  C2.  and  yi  in  Eq.  14  must  be  replaced  by  C/> 

C/+1 ,  and  Yi  ■  Similarly,  1  must  be  replaced  by  i  in  Eq.  17.  The  definition  of  x  and  y  in  terms  of  layer 
1  is  indicated  for  a  profile  consisting  of  a  few  layers.  However,  there  might  be  some  advantage  to 
selecting  a  different  reference  layer  for  complicated  multilayer  profiles.  For  example,  consider  a  major 
refractive  duct  with  a  minor  surface  duct  above.  In  this  case,  we  would  recommend  that  x  and  y  be 
referred  to  the  lower  branch  of  the  major  duct.  With  this  reference,  the  effects  of  the  minor  duct  will 
be  to  change  the  x,y  solution  slightly  from  that  of  the  major  duct  alone.  In  contrast,  the  x,y  solution 
couched  in  terms  of  the  first  layer  will  bear  relatively  little  relation  to  the  solution  for  the  surface  duct 
alone. 

We  next  examine  two  different  implementations  of  the  canonical  eigenvalue  method.  The  simplest 
implementation  is  to  take  a  specific  set  of  profile  parameters,  e.g.,  C,  and  y, .  For  a  general  profile 
consisting  of  n  interfaces  and  m  boundaries,  there  are  2n+m+l  profile  parameters.  We  then  proceed 
to  evaluate  the  coefficients  of  Eq.  36  to  39,  plus  higher  order  terms  necessary,  for  the  given  profile 
parameters. 

This  method  is  simple  but  makes  no  use  of  the  canonical  nature  of  the  method.  The  coefficients 
of  Eq.  36  to  39  depend  only  on  the  parameters  pf  and  S, .  For  a  profile  with  n  interfaces  and  m 
boundaries,  there  are  n  values  of  p,  and  n+m-2  values  of  Si  •  Suppose  now  that  we  specify  a  set  of 
2n+m-2  numerical  values  of  p,  and  Si .  evaluate  the  coefficients  of  Eq.  36  to  39,  and  solve  the 
eigenvalue  equation.  (In  passing,  we  note  that  we  must  specify  p;  rather  than  p?  as  would  be  possible 
in  Eq.  36  to  39  because  p,  appears  in  the  eigenvalue  equation  as  a  coefficient  of  some  of  the  Airy 
functions.)  We  then  have  the  solution  for  arbitrary  values  of  C\>  C2,  and  y,  for  the  profile  class.  The 
canonical  profile  class  is  determined  by  the  set  of  p,  and  S,  and  the  boundary  conditions.  Once  we 
select  a  set  of  values  for  C\,  C2,  and  y, ,  the  remaining  2n+m-2  profile  parameters  for  the  specific 
profile  can  be  generated  from  the  set  of  p,  and  5, . 

This  generation  is  accomplished  by  recursion.  Equation  29  may  be  solved  to  yield  C,t2  as  given  by 

CU2  =  Ci(l(Ci/CUt)2  -  1]  Su j  +  l}-,/2  .  (48) 

We  first  determine  Cj  in  terms  of  the  given  C,.  C2*  and  S2  and  use  Eq.  48  recursively  to  generate  all 
values  of  Cj .  Equation  9  may  be  solved  for  ym  to  yield 

Yu l  =  -Yiopf  ■  (49) 

We  initially  use  Eq.  5  to  generate  y10  in  terms  of  the  given  CP  C2,  and  y1.  We  then  use  Eq.  49 
to  solve  for  y2.  We  then  alternately  use  Eq.  5  and  49  to  generate  yUl  from  the  determined  y,  and 
Cj  and  Ci„  as  previously  established  with  the  use  of  Eq.  48. 

It  is  of  interest  to  compare  the  canonical  eigenvalue  approach  with  the  customary  approach  in 
terms  of  free  variables,  dependent  (constrained)  variables,  and  parameters.  We  will  define  parameters 
as  those  quantities  which  are  fixed  in  the  solution  of  the  eigenvalue  equation.  We  will  define  free 
variables  as  those  that  we  may  choose  independently  of  the  eigenvalue  equation.  Consider  the  general 
profile  of  n  interfaces  and  m  layers.  The  eigenvalue  equation  in  the  customary  approach  contains  the 
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profile  parameters  plus  the  frequency,  for  a  total  of  2n+m+2  parameters  and  one  dependent  variable, 
Cp,  as  constrained  by  the  eigenvalue  equation. 

In  the  canonical  approach,  we  add  the  variables  x  and  y.  Thus,  we  must  account  for  a  total  of 
2n+m+5  values,  which  may  be  either  free  variables,  dependent  variables,  or  parameters.  The  parame¬ 
ters  in  this  approach  are  given  by  the  2n+m-2  values  of  p,-  and  Si-  The  free  variables  are  C\<  C2,  yx , 
and  Jt.  The  dependent  variables  are  y,  /,  and  Cp  as  determined  respectively  by  the  constraints  of  the 
canonical  eigenvalue  equation,  by  Eq.  14,  and  by  Eq.  17.  This  accounts  for  the  requisite  number  of 
n+m+5  variables. 

The  significance  of  the  canonical  eigenvalue  method  is  that  the  number  of  parameters  in  the 
eigenvalue  equation  has  been  reduced  by  four  as  compared  to  the  customary  approach.  The  four 
parameters  that  have  been  removed  are  C\,  Ci>  yx>  and/.  We  indicate  the  mathematical  features  that 
allow  these  to  be  removed.  Equation  29  contains  three  sound  speeds.  This  allows  one  to  treat  the 
canonical  eigenvalue  equation  as  independent  of  Cx  and  Ci-  Here,  Sz  in  the  eigenvalue  equation 
determines  Cj  in  Eq.  48  for  any  choice  of  C\  and  Cz-  Similarly,  y,  is  independent  because  it  is  the 
only  gradient  which  does  not  appear  in  the  definition  of  the  p,  of  the  canonical  eigenvalue  equation. 
The  parameter  /  is  removed  from  the  eigenvalue  equation  by  making  it  satisfy  the  constraint  of  Eq.  14 
in  terms  of  the  removed  parameters  C\<  Cz<  and  yx ,  and  the  independent  variable  x. 

Consider  now  the  one-layer  bounded  profile.  The  four  parameters  appearing  in  the  customary 
eigenvalue  equation  are  reduced  to  zero.  In  this  case,  the  canonical  eigenvalues  are  independent  of  all 
profile  parameters.  The  solutions  apply  to  all  profiles  of  the  class,  where  the  class  is  defined  as  a  one- 
layer  profile  with  specified  free,  rigid,  or  combination  boundary  conditions. 

Consider  next  the  thiee  simple  proai;s  treated  in  Ref.  1.  Two  of  the  three  profiles  were  one-layer 
surface  ducts  bounded  above  by  a  free  or  rigid  boundary.  For  these  profiles,  the  three  parameters  of 
the  customary  approach  were  reduced  to  zero  in  the  canonical  eigenvalue  equation.  For  the 
unbounded  two-layer  refractive  profile,  the  four  parameters  of  the  customary  approach  are  reduced  to 
one  in  the  canonical  approach,  i.e.,  the  parameter  p*.  These  three  simple  profiles  are  special  cases  in 
which  either  x  or  y  alone  but  not  both  appear  in  the  eigenvalue  equation.  Here,  Eq.  14  does  not 
apply  as  the  frequency  is  independent  of  x  or  y  and  the  phase  velocity  for  all  frequencies  may  be  gen¬ 
erated  from  Eq.  17  or  its  y  equivalent. 

The  two-layer  bounded  duct  is  the  simplest  profile  which  exhibits  the  phenomenon  of  double-duct 
propagation.  As  discussed  in  Ref.  1,  we  obtained  numerical  examples  of  this  phenomenon  but  were 
unable  to  cope  with  it  analytically  because  of  the  large  number  of  parameters  in  the  customary 
approach.  The  canonical  eigenvalue  approach  reduces  the  six  parameters  (Cx,  Cj.  C3,  yx ■  yi,  and/) 
for  the  two-layer  bounded  profile  to  two,  viz.,  p;  and  S2.  Our  original  goal  in  the  process  of  develop¬ 
ing  the  canonical  eigenvalue  method  was  to  facilitate  the  theoretical  analysis  of  the  double-duct  propa¬ 
gation.  It  appears  now  that  this  goal  has  been  met,  but  the  degree  of  success  remains  to  be 
demonstrated. 

In  the  case  of  many  profile  layers,  the  advantage  of  the  canonical  forrhulation  in  reducing  the 
number  of  parameters  is  minimal.  For  example,  in  the  case  of  a  four-layer  profile,  the  reduction  from 
10  to  6  parameters  is  not  particularly  helpful.  However,  we  believe  that  the  canonical  formulation  may 
still  have  advantages  over  the  customary  formulation,  even  when  the  solution  involves  the  numerical 
evaluation  of  the  eigenvalue  equation  in  matrix  form.  We  believe  that  the  solution  of  the  eigenvalue 
equation  in  x,  y  space  will  be  more  predictable  than  the  solution  in  Cp  and  /  space. 

We  will  now  determine  if  there  are  any  restrictions  on  the  choice  of  the  set  of  p,  and  Sj  in  the 
canonical  eigenvalue  equations.  We  will  first  do  the  analysis  for  the  two-layer  case. 
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As  will  be  illustrated  in  more  detail  later,  the  eigenvalue  solutions  and  profile  class  are  divided 
into  Case  A  as  given  by  Eq.  16  and  Case  B  as  given  by  Eq.  15.  It  is  convenient  to  treat  these  sepa¬ 
rately.  Consider  first  Case  A,  with  p,  positive.  Here  there  are  two  possibilities.  The  first  is 


C,SC3<Cj.  (50) 

From  Eq.  29,  we  determine  that  when  C3  =  C2,  =  1»  and  when  C3  -  Cx,  S2  =  0.  Thus,  the 

allowable  values  of  Si  under  the  conditions  of  Eq.  50  are 

0SS2<1.  (51) 


The  second  possibility  is 


0  S  Cj  S  C|  <  c2 . 


(52) 


The  counterpart  of  Eq.  52  is 


-  ®  <;  s2  <;  o  .  (53) 

Consider  next  Case  A,  with  pj  negative.  Here  the  only  possibility  is 

C|SC2<Cj<«.  (54) 


The  allowable  values  of  Si  are 


1  <S2<  ®  .  (55) 

Consider  now  Case  B.  The  counterparts  of  Eq.  50,  52,  and  54  are,  respectively 

C2  <  C3  S  Cx  (56) 

Ct<C,SC,<®  (57) 

and 

0  S  C3  S  C2  <  C|  .  (58) 

However,  the  restrictions  on  Sz  are  the  same  as  for  Case  A;  i.e.,  Eq.  51,  53,  and  55,  respectively, 
apply  here. 

We  see  then  that  the  only  restriction  on  p!  and  S2  is  that  S2  must  be  greater  than  1  for  p,<0  and 
less  than  1  for  p2>0.  The  analysis  is  similar  for  p,  and  S,+1.  Here,  the  restriction  is  that  Sui  must  be 
greater  than  1  for  p,<0  and  less  than  1  for  p,>0.  These  represent  very  mild  restrictions.  The  case  of 
p,  negative  is  not  of  great  theoretical  interest  since  it  merely  changes  the  degree  of  slope.  The  case  of 
p,  positive  reverses  the  sign  of  the  slope  and  sets  up  a  double-duct  configuration.  Values  of  5,tI 
between  1  and  -1  will  be  of  most  interest  since  values  of  S,+1  <  -1  may  correspond  to  unrealistically 
low  values  of  sound  speed. 

For  any  values  of  Sf+1,  both  Case  A  and  Case  B  solutions  are  possible.  The  canonical  eigenvalue 
curves  are  divided  into  two  sets— those  for  Case  A  and  Case  B  profiles.  This  works  fine  for  bounded 
profiles.  However,  there  is  a  problem  when  the  deepest  profile  layer  is  unbounded.  As  discussed  in 
Ref.  1,  the  canonical  eigenvalue  approach  was  first  used  on  a  two-layer  refractive  duct  bounded  above 
by  a  free  surface,  with  the  second  layer  unbounded.  The  canonical  eigenvalue  curves  were  valid  for 
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Case  B,  but  not  for  Case  A,  where  the  half  space  has  a  negative  rather  than  positive  gradient.  The 
problem  is  that  the  solution  for  a  positive-gradient  half  space  involves  the  Airy  function  Ai  whereas 
that  of  a  negative  gradient  half  space  involves  the  modified  Hankel  function  h2  ■  Thus,  for  unbounded 
ducts,  the  canonical  eigenvalue  equation  must  take  two  separate  forms— one  for  Case  A  and  one  for 
Case  B. 

For  an  unbounded  duct  with  a  negative  gradient,  the  phase  velocity  has  a  complex  component 
which  results  in  attenuated  modes  associated  with  leakage  out  of  the  channel.  In  such  a  case,  the 
eigenvalues  x  and  y  will  be  complex.  We  will  now  consider  complex  eigenvalues  to  see  what  particular 
problems  are  involved.  From  Eq.  7  and  8  and  the  fact  that  Cp  is  the  only  complex  term  in  the  right 
side,  it  follows  that 


Im  =  Im  fa  ■  (59) 

In  particular,  it  follows  that 

Im  x  =  /my.  (60) 

Equation  60  also  results  in  a  real  frequency  in  Eq.  14. 

We  have  not  proved  that  Eq.  60  holds,  but  we  have  inferred  it  from  the  relationship  between  the 
physical  and  mathematical  variables.  There  must  be  something  about  the  mathematical  form  of  the 
canonical  eigenvalue  equation  that  leads  to  Eq.  60.  This  analysis  is  beyond  the  scope  of  this  report 
and  is  best  left  until  specific  examples  of  complex  eigenvalues  are  considered.  We  can  demonstrate 
that  Eq.  59  holds  if  we  assume  that  Eq.  60  holds.  We  assume  that  Eq.  59  holds  for  /'  and  need  to 
demonstrate  that  it  holds  for  /+1.  From  Eq.  18  and  Eq.  26,  it  follows  that 


Im  o  =  {Mux  +  Nux)  Im  £;0  =  pf  Im  ?/0  .  (61) 


From  Eq.  10,  it  follows  that 


Im  $m.i  =  p/  Im  pto  ■  (62) 

and  we  have  demonstrated  that  Eq.  59  holds  for  j+1.  In  passing,  it  is  of  interest  to  note  that  the 
imaginary  component  of  the  Airy  function  arguments  only  involves  Im  x  and  the  set  of  pj . 

Other  than  the  derivation  of  Eq.  60,  we  see  no  particular  obstacles  to  the  application  of  the 
canonical  approach  for  complex  eigenvalues.  The  eigenvalue  curves  will  consist  of  plots  of  /?1  y  vs  XI  x 
and  Im  y  vs  /?1  x.  These  should  suffice  when  Im  y  is  small  relative  to  /?!  y,  which  is  the  usual  case. 

We  have  just  discussed  a  case  where  the  profile  parameters  are  real  but  the  eigenvalues  of  both 
the  customary  and  canonical  approach  are  complex.  We  now  point  out  an  application  where  the  pro¬ 
file  parameters  are  complex.  As  discussed  in  Ref.  3,  one  of  the  techniques  for  modeling  bottom  loss  is 
to  introduce  a  sedimentary  bottom  layer  with  complex  values  of  C;  and  y,  in  Eq.  1.  Here,  complex 
parameters  are  introduced  into  the  customary  eigenvalue  equation,  with  the  result  that  Cp  has  an 
additional  complex  component  due  to  attenuation  in  the  bottom  sediment. 

The  canonical  approach  can  simplify  this  analysis  if  we  define  x  and  y  as  associated  with  the  sedi¬ 
mentary  bottom  layer.  We  thus  eliminate  the  effect  of  the  complex  parameters  on  the  canonical  eigen¬ 
value  equation.  The  effect  of  the  complex  parameters  is  contained  in  the  equivalent  of  Eq.  17  for  the 
sedimentary  layer.  One  can  then  make  a  parametric  study  of  how  complex  C/  and  y,  for  the  sedi¬ 
mentary  layer  affect  the  model  attenuation  for  the  profile  class.  In  the  case  of  a  bounded  profile,  x  in 
Eq.  17  can  be  real.  In  the  case  of  an  unbounded  profile,  x  will  have  an  imaginary  argument  that  also 
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contributes  to  the  attenuation.  In  the  canonical  approach,  the  attenuation  associated  with  complex  x  is 
clearly  separated  from  that  caused  by  complex  profile  parameters.  By  contrast,  in  the  customary 
eigenvalue  approach,  the  two  attenuation  factors  are  scrambled  together,  yielding  only  the  combined 
complex  component  of  Cp- 

There  are  several  details  to  be  worked  out  before  this  application  can  be  implemented.  We  would 
like  the  values  of  p,  and  S,  to  be  real,  although  it  is  not  entirely  clear  that  this  can  be  done  when 
some  of  the  profile  parameters  are  complex.  We  also  need  to  treat  the  case  of  a  sound-speed  discon¬ 
tinuity  at  the  water-sediment  interface.  This  will  introduce  some  complication,  which  needs  to  be 
solved.  Finally,  the  presence  of  complex  parameters  in  Eq.  14  will  result  in  an  imaginary  component 
in  the  frequency  which  needs  to  be  interpreted. 

GROUP  VELOCITY  AND  EIGENFUNCTIONS 

In  this  subsection,  we  investigate  group  velocity  and  eigenfunctions  to  determine  how  they  can  be 
formulated  in  the  canonical  approach. 

Consider  first  the  group  velocity.  Equation  33  of  Ref.  1  gives  an  expression  for  the  group  velocity 
in  terms  of  the  phase  velocity  as 

Cg  =  Cp  [l-f(dCp/df)/Cp]-'  .  (63) 

We  now  develop  an  expression  for  Cg,  using  the  canonical  approach.  Our  first  step  is  to  square  Eq. 
17,  transpose  the  bracketed  term  to  the  left  side,  differentiate  with  respect  to  /,  and  simplify  to  obtain 

2 (dCp/df)  C]C’PX  +  2Cpf  _1  [1  -CjCp)1]  3'1  (64) 

-  Cjf  '2/ V2'3  |y,|2/3  (dx/df)  =  0  . 

It  follows  that 


- f {dCp/df) C~p  =  -(Cp/C,)2  [1-(C,/CP)2]3-1  (65) 

+  (Cp/C,)2  2 ->/ VV2/*  Itt | 2/3  (dx/df)  . 

We  next  take  the  2/3  power  of  Eq.  14  and  differentiate  with  respect  to  /  to  obtain 

2/-./33-i  =  -L-'  (dx/df)(C\C?-\yx  |y,|2/V2/3  (66) 

where 

L  =  (l-dy/dx)-'  .  (67) 

We  solve  Eq.  66  for  dx/df,  substitute  into  Eq.  65,  and  simplify  to  obtain 

-f(dCp/df)C-p'  =  -  (Cp/Ci)2  (1  -  (Ci/Cp)2)  3'1  +  (Cp/C,)2  3-'L(C?C-22-l)  (68) 

Equation  63  may  now  be  written  as 

Q  =  Cp  1 1  +  (Cp  -  C2)/3C2  +  Cp(Cj2  -  Cj2)L/3]'1  .  (69) 

The  application  of  Eq.  69  starts  with  the  evaluation  of  L,  which  we  refer  to  as  the  canonical 
group  velocity  factor.  When  the  eigenvalue  curves  are  generated,  we  also  generate  plots  of  L  vs  x. 
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Equation  14  then  gives  the  frequency  for  x  and  the  desired  profile  parameters.  Equation  17  generates 
the  phase  velocity  and  we  can  then  evaluate  Eq.  69,  plotting  it  as  a  function  of  frequency. 

One  cannot  appreciate  the  elegance  of  Eq.  69  unless  one  has  evaluated  group  velocity  for  the 
conventional  approach.  As  discussed  in  Ref.  4,  this  can  be  done  by  involved  matrix  methods.  For 
profiles  with  many  layers,  dyldx  will  need  to  be  done  by  a  similar  matrix  method.  However,  the  for¬ 
mulation  should  be  somewhat  easier  to  implement  than  the  conventional  approach.  In  the  case  of  the 
simple  ducts  of  Ref.  1,  where  y  or  x  is  missing,  the  last  term  of  Eq.  69  is  suppressed.  Reference  1 
demonstrates  that  Eq.  69  is  the  same  result  that  we  would  get  from  the  ray  theory  result 

Cg  =  R  /T  (70) 

where  R  and  T  are  the  cycle  range  and  travel  time,  respectively. 

In  the  case  of  eigenfunctions,  the  same  general  procedures  apply  here  as  in  the  conventional 
approach.  Here,  the  depth  dependence  of  Eq.  3  may  be  written  as 

(Z)  =  Sit  +  ai(Z-Zj)  .  (71) 

Not  very  much  simplification  here.  The  depends  on  the  profile  class.  However,  a/  depends  on  the 
particular  profile. 

It  is  now  convenient  to  introduce  some  new  notation.  We  define  the  following: 


Fi(-S„)  =  DjAi(-  Sii)  +  EiBK-Zn)  .  (72) 

F«(-$«o)  =  DiAi(-£i 0)  +  EiBK-b o)  ,  (73) 

Fi'  =  dF/dt  ,  (74) 

Fi'(-Sn)  =  -ZW(-S/,)  -  £,*/'(-&,).  (75) 

and 

Fi'l-ti o)  =  -fW(-Sio)  -  EiBi'(-Zio)  •  (76) 

We  note  that  the  continuity  conditions  at  the  interface  result  in 

Fht  (-&♦,.,)  =  Fj(-S;o)  (77) 

and 

F',>, (-£/♦,.,)  =  -p/F,'(-S,o)  .  (78) 

Our  next  step  is  to  evaluate  the  coefficients  £>,  and  E\  of  Eq.  2.  Let  us  assume  that  we  know  £>, 
and  Ei  and  wish  to  evaluate  D,»,  and  £,■*, .  Equations  77  and  78  represent  a  pair  of  linear  equations 
in  Djtl  and  Ejtl  in  which  the  right  side  is  known.  The  solution  is 

Du i  =  *  [F/(-S/0)S<'(-pjS,o)  -  PiFi'(-S,o)Fi(-p?5,o)l  (79) 

and 

Eui  =  -»[F<(-£/oM<'(~P??/o)  -  P/F,’(-UM/(- pj&o)]  •  (80) 
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The  tt  factor  in  Eq.  79  and  80  is  the  reciprocal  of  the  Wronskian  of  Airy’s  differential  equation. 

We  start  the  process  at  the  ocean  surface.  For  a  free  surface 

£>i  =  fli'(-Sn)  and  Ei  =  -Ai(-?„)  (81) 

satisfies  the  boundary  condition.  The  counterpart  of  Eq.  81  for  a  rigid  surface  is 

Dy  =  Bi'(-t ,,)  and  E\  =  .  (82) 

We  initially  use  Eq.  81  or  82  to  evaluate  £i(-?io)  and  Fx'{-  £10)-  We  then  use  Eq.  79  and  80  to 
evaluate  Z>2  and  £2.  We  then  use  these  to  evaluate  £2(-£20)  and  £2'(-£2o)  and  use  Eq.  79  and  80  to 
evaluate  D3  and  £3.  This  process  is  repeated  until  Dj  and  Et  are  evaluated  for  the  last  profile  layer. 

In  passing,  we  note  that  the  above  procedure  satisfies  all  boundary  and  interface  conditions  except 
the  last  one.  In  the  case  of  a  boundary  at  the  bottom  of  layer  /,  this  procedure  does  not  satisfy  the 
bottom  boundary  condition,  because  Dj  and  Ej  are  already  determined.  Similarly,  in  the  case  of  an 
unbounded  layer  with  positive  gradient  at  the  last  interface,  one  condition  will  not  be  satisfied.  What 
happens  is  that  the  eigenvalue  equation  itself  guarantees  the  satisfaction  of  this  condition. 

The  procedure  outlined  in  Eq.  77  to  82  is  essentially  the  same  as  for  the  customary  approach. 
There  is  one  interesting  feature.  The  evaluation  only  involves  quantities  which  are  available  from  the 
canonical  eigenvalue  solution.  What  this  means  is  that  the  coefficients  D,  and  £,  can  be  specified  for 
the  entire  profile  class. 

We  now  turn  to  the  problem  of  normalization.  The  eigenvalue  equation  ensures  a  solution  to  the 
homogeneous  equations  which  satisfies  the  boundary  and  interface  conditions.  However,  this  solution 
is  not  unique,  since  any  constant  times  the  solution  is  also  a  solution.  A  unique  multiplicative  constant 
is  determined  by  a  normalization  procedure.  We  define  the  normalization  factor  Dn  as 


°"  =  Z  7  F<(Z)  dZ  •  <83> 

/7l  JZ> 

The  normalized-mode  depth  functions  are  given  by 

U„(Z0)  U„(Z)  =  F„(Z 0)  F„{Z)/D„  .  (84) 

Here  n  is  the  mode  number,  Z0  is  the  source  depth,  and  Z  is  the  receiver  depth.  The  F„(Z0)  or 
£„(Z)  °f  Eq-  is  the  £j  of  Eq.  2  for  the  layer  in  which  the  source  or  receiver  is  located,  as  evalu¬ 
ated  for  mode  n. 

We  start  the  solution  Eq.  83  by  noting  that  Eq.  2  is  a  solution  of  Airy’s  differential  equation  as 
given  by 

dlF/d?  -(-£)£  =  0  .  (85) 

We  next  note  that 

F2(-0  =  d[£F2  +  (dF/d02]/d£  ■  (86) 
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We  obtain  F2  if  we  carry  out  the  derivative  in  the  right  side  of  Eq.  86  and  make  use  of  Eq.  85.  It 
follows  then  that 


r  Zi+1  Z/+  l 

F](Z)dZ  =  aj'[ZF2  +  (dFi/dt)2)  .  (87) 

J  Zi  Zi 

Equation  87  makes  use  of  the  relationship  obtained  from  Eq.  3  of 

dt/dZ  =  at  .  (88) 

Our  goal  here  is  not  only  to  solve  for  Dn  but  to  cast  it  into  the  form 

Dn  *  a T1  De  (89) 

where  a;  is  given  by  Eq.  4  and  Dc  depends  only  on  quantities,  which  are  available  from  the  canoni¬ 
cal  eigenvalue  equation.  We  will  refer  to  Dc  as  the  canonical  normalization  coefficient.  The  advan¬ 
tages  of  this  particular  breakdown  will  soon  be  apparent. 

For  purposes  of  analysis,  it  is  convenient  to  express  Dc  as 

/- i 

Dc  =  Ds  +  Db  +  £  Dm  ■  (90) 

im  1 

The  quantities  Ds  and  £>fc  represent  the  contributions  at  the  surface  and  bottom  boundaries.  The 
quantity  D„/  represents  the  contribution  to  Db  at  interface  i.  This  contribution  is  summed  over  all 
interfaces,  which  total  7-1  for  a  bounded  profile  with  7  layers. 

From  Eq.  11,  it  follows  that 


oil  i  =  -Pioj1  ■  (91) 

Hence 

oj'  =  Dpi  a-'  (92) 

where 

<-i 

Dpi  =  ft  ("P7>  ’Dp  1  =  1-  (93) 

7-i 

Equation  92  follows  from  the  repeated  use  of  Eq.  91  and  expresses  ajl  as  aj1  times  the  Dpi  of  Eq. 
93,  which  is  a  product  function  of  the  negative  of  py  for  all  /  less  than  /.  From  Eq.  87,  we  see  that 
aj1  is  a  factor  associated  with  the  layer  7  contribution  to  D„ .  Equation  92  allows  us  to  express  Dn  in 
the  form  of  Eq.  89. 

We  now  evaluate  the  contribution  to  Dc  at  the  surface  boundary.  In  the  case  of  a  free  surface, 
Ft  =  0  and  the  lower  limit  of  Eq.  87  leads  to 

D,  =  —  ji~2  (94) 
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where  is  the  Wronskian  for  Airy’s  differential  equation.  In  the  case  of  a  rigid  surface,  F\  =  0, 
and  the  lower  limit  of  Eq.  87  can  be  simplified  to 

Ds  =  -?n*-2  •  (95) 

We  next  proceed  to  evaluate  the  contribution  at  interface  i+1.  Equation  87  is  by  layers,  not  inter¬ 
faces.  At  interface  i+1,  we  have  the  upper  limit  contribution  of  Eq.  87  for  layer  <  and  the  lower  limit 
contribution  of  Eq.  87  for  layer  i+1.  We  make  use  of  Eq.  77  and  78  to  express  the  layer  it  1  contri¬ 
bution  in  terms  of  the  layer  i  contribution.  The  application  of  Eq.  87  then  leads  to 

aj1  Dni  =  (a]% 0  -aj*t  5m .,)F?(-fco)  +  (aj1  -  aJ^F'fC-U)  •  (96) 

Equation  10  and  11  allow  a  further  simplification  to 


Dni  =  Dpi  (1  +p?)[5/o^?(-5/o)  +  F']  (-Ml  ■  (97) 

In  the  case  of  an  unbounded  medium,  the  contributions  to  Eq.  87  at  infinity  are  taken  as  zero. 
Thus,  there  is  nothing  more  to  add  and  Db  of  Eq.  90  is  zero.  In  the  case  of  a  medium  with  lower 
bound,  we  must  evaluate  the  contribution  at  this  boundary.  In  the  case  of  a  free  boundary  at  the  bot¬ 
tom  of  layer  /,  we  obtain 


Db  =  Dpl  F')&, 0)  .  (98) 

In  the  case  of  a  rigid  bottom,  we  obtain 

Db  =  Dp/5/0Fj(-5/o)  •  (99) 

All  components  in  Eq.  94,  95,  and  97  to  99  are  expressed  in  terms  of  quantities  that  are  avail¬ 
able  from  the  canonical  eigenvalues  and  canonical  profile  parameters.  The  particular  profile  parame¬ 
ters  y,  and  C\  and  the  frequency  enter  into  D„  through  the  simple  multiplicative  factor  aj1  in  Eq. 
89.  The  bulk  of  the  normalization  evaluation,  Dc,  can  be  evaluated  for  the  canonical  profile  class. 
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SECTION  2.  OUTLINE  OF  EVALUATION  BY  THE 
CANONICAL  APPROACH 

This  section  outlines  the  general  steps  necessary  for  the  evaluation  of  normal-mode  fields  by  the 
canonical  approach.  It  also  serves  as  a  summary  of  the  method. 

We  start  with  a  bounded  sound-speed  profile  class  of  /  layers.  The  class  is  characterized  by  (/-l) 
values  of  p,  and  £,  and  by  a  given  pair  of  boundary  conditions.  The  I  X  I  eigenvalue  matrix  of  the 
customary  approach  is  formed  by  using  the  boundary  and  interface  conditions.  The  arguments  of  the 
Airy  functions  in  layer  i  are  given  genetically  by  -  £/i  at  the  upper  interface  and  -  £/o  at  the  lower 
interface. 

The  next  step  is  to  use  Eq.  32  and  33  to  replace  these  arguments  by  linear  functions  of  x  and  y. 
The  coefficients  of  these  linear  functions  are  functions  p,-  and  Si  and  are  generated  by  recursion  with 
Eq.  18.  Equations  34  and  39  give  these  coefficients  for  the  first  three  layers.  At  the  end  of  this  proc¬ 
ess,  we  have  the  canonical  eigenvalue  equation,  which  we  may  describe  generically  as 

G(x,  y,  pi,  Si)  =  0  .  (100) 

The  solution  of  this  equation  gives  the  canonical  eigenvalues,  y  as  a  function  of  x.  These  eigenvalues 
apply  to  the  canonical  class  of  the  sound-speed  profile  defined  by  the  set  of  p,  and  5,  and  the  given 
boundary  conditions. 

The  next  stage  in  the  process  is  to  generate  the  eigenvalue  curves,  y  as  a  function  of  x,  for  the 
given  set  of  p,-  and  Si  ■  This  is  accomplished  by  iteration  using  Newton’s  method  where  dG/dx  or 
dG/dy  are  the  required  derivatives.  These  partial  derivatives  are  also  used  to  obtain  dy/dx  so  that  the 
canonical  group  velocity  factor.  L  of  Eq.  67,  can  be  evaluated  for  later  use.  The  canonical  frequency 
function,  (y~x)3/2,  and  the  canonical  normalization  coefficient,  Dc  of  Eq.  90,  are  also  evaluated  for 
later  use. 

All  the  quantities  generated  thus  far  apply  to  the  canonical  profile  class.  To  proceed  further,  we 
must  select  the  parameters  Cp  C2.  and  yi  of  the  first  layer  of  the  particular  profile  we  wish  to  evalu¬ 
ate.  These  parameters  plus  the  set  of  p,-  and  Si  completely  define  the  sound-speed  profile. 

The  next  step  is  to  use  Eq.  14,  the  canonical  frequency  function,  Cp  C2.  and  yi  to  determine 
the  frequency,  /,  as  a  function  of  x.  The  next  step  is  to  use  Eq.  17,  x,  /,  Cp  and  yi  to  determine  the 
phase  velocity,  Cp,  as  a  function  of  /.  The  next  step  is  to  use  Eq.  69,  L,  Cp<  Cp  and  C2  to  deter¬ 
mine  the  group  velocity,  Cg>  as  a  function  of  /.  The  final  step  is  to  use  Eq.  89,  £>c.  /,  Cp  and  y,  to 
determine  the  eigenfunction  normalization  coefficient,  Dn,  as  a  function  of  frequency.  The  latter  step 
does  not  depend  on  the  phase  or  group  velocity  and  could  be  carried  out,  if  so  desired,  before  Cp 
and  Cg  are  evaluated. 

The  depth  functions  of  Eq.  84  may  now  be  evaluated  as  a  function  of  depth  to  give  the  standing 
waves.  The  range  dependence  is  given  by  H\{wr/Cp)t  where  H\  is  a  Hankel  function  of  the  second 
kind.  The  acoustic  field  is  then  evaluated  by  multiplying  the  depth  functions  by  the  range  dependence 
and  summing  over  the  modes. 

Reference  2  applies  the  various  stages  of  this  approach  to  the  class  of  one-layer  bounded  profiles. 
Three  pairs  of  boundary  conditions  are  treated— a  free  surface  and  rigid  bottom,  a  free  surface  and 
free  bottom,  and  a  rigid  surface  and  rigid  bottom.  Numerical  results  are  presented  to  illustrate  the 
desired  outputs  and  other  key  elements  of  the  evaluation. 
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